GENERATION OF PULSAR GLITCHES: A SUPERFLUID 

CORE MODEL a 



A. SEDRAKIAN & J. CORDES 

Center for Radiophysics and Space Research, Cornell University, 
Ithaca, NY 1+853 

We show that the neutron star's crust-core interface acts as a potential barrier on 
the peripheral neutron vortices approaching this interface and thus prevents their 
continuous decay on it in the course of equilibrium state deceleration. The barrier 
arises due to the interaction of vortex magnetic flux with the Meissner currents set 
up by the crustal magnetic field at the interface. When the non-balanced part of 
the Magnus force reaches the value at which the vortices are able to annihilate at 
the interface, the rapid transfer of angular momentum from the superfluid spins-up 
the observable crust on short dynamical coupling times. 



Models of generation of pulsar glitches are supposed to explain the fol- 
lowing observational facts: (i) the short spin-up time scales, which are less 
than 120 s in the Vela pulsar 0833-45 and less than an hour in the Crab pul- 
sar 0531+21; (ii) the magnitudes of the jumps in the rotation and spin down 
rates, Au/v ~ 10~ 8 — 10~ 6 and Az>/z> ~ 10~ 3 — 10~ 2 , respectively; and (hi) the 
origin of the instability driving a glitch along with characteristic intervals be- 
tween glitches (typically of the order of several months to years) n A number of 
existing generic models invoke as trigger crust- and core-quakes a (discontinu- 
ous adjustments of the solid crust to the gradually changing oblateness of the 
star as it spins down) ; spontaneous quantum transitions of rotating superfluid 
between quasistationary mgenstates corresponding to different eigenvalues of 
total angular momentumB; collective unpinning of a large number (~ 10 13 ) of 
vortices in the neutron star crustsBu and thermal instabilities^. The increasing 
bulk of observational evidence provides good chances of discrimination between 
the theoretical models in the future. 

Here we shall give a brief account of a new trigger mechanism for generation 
of pulsar glitches in neutron star's superfluid core. The complete discussion will 
be given in ref. 6. In the interjump epoch a neutron star is decelerating; conse- 
quently the vortex lattice in the superfluid core is expanding and the peripheral 
vortices attempt to cross the crust-core boundary. The crust-superfluid core 
interface acts as a potential barrier on the vortices in the superfluid core that 
approach this boundary. The barrier arises due to the magnetic interaction 
between the crustal magnetic field Hq (which penetrates the superconducting 
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core exponentially within a scale S p - the penetration depth) and quantum 
vortices whose magnetic field is governed by the generalized London equation 



5~ 2 V x (V x B v ) + B v = u p $ d<- 2 \r- r p ) + v n <P^ 2 \r- r n ). 



(1) 



Here z?'s are circulation unit vectors; subscripts p and n refer to protons and 
neutrons; $o is the flux quantum carried by proton vortices; and $1 is the 
non-quantized flux of neutron vortices due to the entrainment effect (i.e. the 
effect of superconducting proton mass transport by the neutron superfluid 
circulation, see ref. 7). The total field B is the superposition of B v of eq. ([[]) 
and crustal filed B cr = _ff exp(— f •n/S p ), (ft being the normal of the interface). 

Suppose that the interface is the (yz) -plane of a Cartesian system of co- 
ordinates and the vectors of vortex circulation are in the positive z-direction. 
The half-plane x < corresponds to the crust while x > corresponds to the 
superfluid core. The knowledge of magnetic field distribution with the bound- 
ary condition at the interface B z = H , allows one to calculate the relevant 
part of the Gibbs free energy of the system G = F — (47r) _1 J B ■ H dV, where 
the free-energy is 
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The force per single vortex of effective flux $*(= $0; $1) derived from G is 
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Here K\ is the modified Bessel function. The first term in equation (|^) is 
the repulsive force acting between the vortex magnetic flux and the crustal 
magnetic field. It can also be interpreted as a Lorentz force resulting from 
superposition of velocity fields of the vortex and the surface Meissner currents. 
The second term is the attractive force acting between the vortex and the 
interface. For large distances the repulsive term dominates; (the attractive 
one goes to zero faster since K\(2x) cx yJn/Ax e~ 2x for x — » 00). For small 
distances the second (attractive) term in eq. ([|) dominates. Thus the repulsive 
part of the vortex - crust-core interface interaction, which dominates at large 
distances, acts as a potential barrier on a vortex approaching the boundary. 
The Magnus force is then balanced by the vortex-interface force. When the 
disbalance between this forces reaches the value at which the vortices annihilate 
at the interface, the resultant rapid transfer of angular momentum from the 
superfluid spins- up the crust. 
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Let us next estimate the interaction. For relevant densities, p ~ 2 x 10 
g cm~ 3 , we have S p ~ 100 fm (e. g. ref. 7) and assuming a conventional 
value for the crustal magnetic field Hq = 10 12 G, we find a maximal repulsive 
force _f max = 3.13 x 10 12 dyn cm -1 . In general, the magnitude of the crustal 
magnetic field at the crust-core interface for different objects can vary in a 
reasonable range 10 9 < Ho < 10 13 G. We find that, the maximal force increases 
(decreases) by two orders of magnitude when the crustal magnetic field is 
increased (decreased) by an order of magnitude. 

Further progress needs to specify the ground state structure of vortices 
in the superfluid core. For the vortex cluster model of ref. 7, 8 (no resid- 
ual field when the superconducting state sets on or a complete Meissner ex- 
pulsion of residual field) the number of proton vortices per neutron vortex 
interacting with the interface is N < 283 and the maximal force is /™ ax — 
N / max = 8.9 x 10 14 dyn cm -1 . It should sustain the Magnus force excess 
(i.e. the part which is non-balanced along the normal n) acting on a neutron 
vortex Sf M = 3.23 x 10 17 (<5w s /s _1 ) dyn cm -1 , where 8u> s is the angular ve- 
locity difference between the superfluid and the normal components.^] From 
the balance condition /™ ax = Sf M , the maximal departure that can be sus- 
tained by the boundary force on the cluster is <5w™ ax ~ 0.003 s _1 . Then the 
angular momentum conservation in a Vela-type glitch implies that the ratio 
of the moment of inertia of the superfluid region to the normal component 
should be I s /I n = 0.023. This value is close to a previous estimate of the 
ratio I s /I n = 0.02 for a superfluid shell at the crust-core boundary with short 
dynamical coupling times (< 120 s)u. 
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b Thc neutron star model used in our estimates is discussed in detail in ref. 8. 
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